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A supersymmetric one-dimensional matrix procedure similar to relationships of the 
same type between Dirac and Schrodinger equations in particle physics is described at 
the general level. By this means we are able to introduce a nonhermitic Hamiltonian 
having the imaginary part proportional to the solution of a Riccati equation of the Wit- 
ten type. The procedure is applied to the exactly solvable Morse potential introducing 
in this way the corresponding nonhermitic Morse problem. A possible application is 
to molecular diffraction in evanescent waves over nanostructured surfaces. 
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Un procedimiento matricial uni-dimensional supersimetrico similar a relaciones del 
mismo tipo entre las ecuaciones de Dirac y Schrodinger que usamos recientemente 
para el oscilador armonico clasico es presentado de manera concisa en terminos gen- 
erales. El aspecto nuevo es el uso de parametros constantes por medio de los cuales el 
Hamiltoniano se vuelve no hermftico con parte imaginaria proporcional a la solucion 
de una ecuacion de Riccati de tipo Witten, que es caracterfstica para el metodo su- 
persimetrico. El factor de proporcionalidad contiene los parametros mencionados. 
Aplicamos esta tecnica algebraica al oscilador cuantico no armonico de Morse obte- 
niendo una forma no hermi'tica. Una posible aplicacion es a la difraccion molecular en 
ondas evanescentes sobre superficies nanoestructuradas. 
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1 Introduction 



We have recently elaborated on an interesting way of introducing imaginary parts (non- 
hermiticities) in second order differential equations starting from a Dirac-like matrix 
equation QHJ. The procedure is a complex extension of the known supersymmetric 
connection between the Dirac matrix equation and the Schrodinger equation A 
detailed discussion of the Dirac equation in the supersymmetric approach has been 
provided by Cooper el al in 1988, who showed that the Dirac equation with a Lorentz 
scalar potential is associated with a susy pair of Schrodinger Hamiltonians. In the 
supersymmetric approach one uses the fact that the Dirac potential, that we denote 
by R, is the solution of a Riccati equation with the free term related to the potential 
function U in the second order linear differential equations of the Schrodinger type. 
Indeed, writing the one-dimensional Dirac equation in the form 

[ap + Pm + I3R{x)\iJj(x) = Etp(x) (1) 

where c = % = 1, p = ~id/dx, m (> 0) is the fermion mass, and R(x) is a Lorentz 
scalar function representing the potential in which the relativistic particle moves. The 

wavefunction is a two-component spinor ( ^ ) and the a and (3 matrices are the 

following Pauli matrices 

°v = ( J 0* ) and ax= ( 1 ) ' 
respectively. Writing the matrix Dirac equation in coupled system form leads to 

[D x + m + R]ipi = E%I> 2 (2) 

[-D x + m + R]ip 2 = E^! . (3) 
By decoupling, one gets two Schrodinger equations for each spinor component, respec- 
tively 

Hixjji = [ — D 2 + ?7jj tpi — eipi , e = E 2 -m 2 , (4) 
where the subindex i = 1,2, and 

Ui(x) = ((m + R) 2 - m 2 =f dR/dx) . 
One can also write factorizing operators for Eqs. 

A ± ^±D x + m + R (5) 

such that 

Hi = A~A + - m 2 , H 2 = A + A~ - m 2 . (6) 
However, we have employed a so-called complex extension of the method by which 
we mean that we considered the Dirac potential R as a purely imaginary quantity 
implying that the Schrodinger potentials Ui are complex, and, as such, we deal with 
nonhermitic problems. We considered previously the cases of the classical harmonic 
oscillator and Friedmann-Robertson- Walker barotropic cosmologies, which correspond 
to the very specific situation in which the Dirac mass parameter that we denoted by 
K was treated as a free parameter equal to the Dirac eigenvalue parameter E. This 
is equivalent to Schrodinger equations at zero energy, e = 0. On the other hand, it 
is interesting to see how the method works for negative energies, i.e., for a bound 
spectrum in quantum mechanics. In this paper, we first briefly describe the method 
and next apply it to the case of Morse potential obtaining a nonhermitic version of 
this exactly-solvable quantum problem. 
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2 Complex extension with a single K parameter 

We consider the slightly different Dirac-like equation with respect to Eq. Q 

■D K W=[a v D x + a w {iR + K)]W = KW , (7) 

where K is a (not necessarily positive) real constant. In the left hand side of the 
equation, K stands as a mass parameter of the Dirac spinor, whereas on the right 
hand side it corresponds to the energy parameter. R is an arbitrary solution of the 
Riccati equation of the Witten type |lj 



R' ± R 2 



u , 



(8) 



of mass K at the fixed energy E = K but in a purely imaginary 



where u is the real part of the nonhermitic potential in the Schrodinger equations we 
get. Thus, we have an equation equivalent to a Dirac equation for a spinor W = 
!>i A _ / w f 
!>o ) ~ \ w b 

potential (optical lattices). This equation can be written as the following system of 
coupled equations 

i£> B< £i + (iR + K)<l>i = K(f>o (9) 

— iD x <f>o + (iR + K)4>o = K0i . (10) 

The decoupling of these two equations can be achieved by applying the operator 
in Eq. 1101 to Eq. © . For the fermionic spinor component one gets 



Dl<t> x - R 2 - D x R-i2KR 
whereas the bosonic component fulfills 



D 2 <f>o 



R" 



D x R-i2KR 











(11) 



(12) 



This is a very simple mathematical scheme for introducing a special type of nonher- 
miticity directly proportional to the Riccati solution. The factorization operators can 
be written in this case in the form 



.4 d 



±iD x + K + iR 



(13) 



that allows to write the fermionic equation 1111 as H\4>\ = (A A + — K 2 )<f)\ = and 
the bosonic equation 1121 as H2(j)o = (A + A~ — K 2 )<f>a — 0. 



3 Complex extension with parameters K and K'. 

A more general case in this scheme is to consider the following matrix Dirac-like 
equation 



-i 

1 



Ox 



1 

1 



iR + K 
iR + K 



wi 

W 2 





°) 


( ^ ) 


V o 


K j 


\ w 2 / 



(14) 
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The system of coupled first-order differential equations will be now 

- iD x + iR + K w 2 = K'wi 
iD x + iR + K wi = k'w 2 
and the equivalent second-order differential equations 



Ac Wi + 



± D X R + 2iKR + (K 2 - K' 2 ) - R 2 



Wi — 



(15) 
(16) 

(17) 



where the subindex i = 1, 2 refers to the fermionic and bosonic components, respec- 
tively. 

Again, introducing the same factorization operators as for the single K case, i.e., 
A ± = ±iD x + K + iR, one can write Eq. 1171 in the Schrodinger-like form 



fliwi = {A~A + - K 2 )w 1 = (K' 2 - K 2 )wi 



(18) 



and 



H 2 w 2 = (A + A~ - K 2 )w 2 = (K 2 - K 2 )w 2 , (19) 

for the fermionic and bosonic components, respectively. These forms are useful for 
quantum mechanical applications; see the next section for one of them. 



4 Application to the Morse potential 

This potential is frequently used in molecular physics in connection with the dissoci- 
ation and vibrational spectra of diatomic molecules. In this case, the Riccati solution 
is of the type 

R{x) = A - Be~ ax , (20) 
Therefore, the second-order fermionic differential equation will be 

D 2 wi + [- (Be -201 - C ie - ax ) + {K 2 - K' 2 ) ~ A 2 

+ 2iK(A - Be~ ax )] wi = (21) 

where B = B 2 , and Ci = B(2A + a). 

The solution is expressed as a superposition of Whittaker functions 

ax/2 x r ( 2B —ax\ , rt ax/ 2 frr / 2B —ax\ 

wi = aie ' M K1 , M ( — e I +/?ie ' W K1 , fl i—e I (22) 

- 1 = ^(2+i"^)and^4(4^- l ^) 1/2 . 
The bosonic equation reads 

D 2 w 2 + [- (Be~ 2ax - C 2 e' ax ) + (K 2 - K' 2 ) - A 2 

+ 2iK(A - Be~ ax )] w 2 = (23) 

where B = B 2 , and C 2 = B(2A - a). 

The solution is a superposition of the following Whittaker functions 

w 2 = a 2 e ax/2 M K2 ^ e-") + f3 2 e ax/2 W K2 ,^ e" a *) (24) 
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where K2 = ^ (2 — — an d the subindex is unchanged. 

If we now place ourselves within the quantum mechanical (hermitic) Morse problem 
we should take 02 = and K — in order to achieve the exact correspondence with the 
bound spectrum problem and eliminate the non-hermiticity. Moreover, the following 
well-known connection with the associated Laguerre polynomials 

M f+n+iif ( 2/ )= y 4 i e -^(y) , V = ^r»" (25) 

can be used in our case with the following identifications 

p K' T ., , . . 

I = — , K ={A — an) 

Then we can write the solution of the hermitic bosonic problem in the well-known 
form 

W2Ay) = Q2 (M) * y^e-^L^- n \y) . (26) 

If we want to approach the nonhermitic problem we define by analogy with Eq. 125H 

M K2lfI (y)=y^e-^ 2 L^ , (y) , (27) 

where K2 and fx are the complex parameters mentioned before and the symbol cor- 
responding to the associated Laguerre polynomial representing now a Laguerre-like 
function introduced by definition through Eq. 12711 . The wavefunction of the nonher- 
mitic problem can be written as follows 

W2,„ „fcerm(?/) = Oi 2 ( — ) 6~ V ' /2 ! (ij) ■ (28) 

In the case of the nonhermitic fermionic problem, the formulas are similar with 
the replacement of ki by k\. Thus: 

wi,, lon h erm (y) = qi ( — \ y^e~ vl2 L 2>1 , (y) . (29) 
V a / Ki-M— 2 

In conclusion, the supersymmetric connection between Dirac-like equations and 
Schrodinger equations, in a simple complex extension form, has been applied here in 
the quantum context of the Morse potential. However, in the pure quantum case, 
the results of this note seem to be only of mathematical interest. A natural question 
is how different types of nonhermiticities can be engineered. As we noticed in our 
previous research |T], physical optics is closer to real applications. In particular, one 
can think to the diffraction of diatomic molecules in evanescent fields because such 
fields have imaginary wavenumbers and we know that Schrodinger equations are sim- 
ilar to Helmholtz equation in the paraxial approximation. A specific experimental 
setup could be very similar to that discussed recently by Leveque and collaborators 
in their study of diffractive scattering of cold atoms from an evanescent field, spa- 
tially modulated by an array of nanometric objects with high index of refraction and 
subwavelength periodicity deposited on a glass surface. The evanescent wavefield is 
created by a totally internally reflected laser beam and is strongly modulated by the 
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configuration of the nanostructure. The calculations are not easy as one should tackle 
Hclmholtz equations in complicated geometries. The task is to obtain the configura- 
tion of the nanostructure that is able to produce the evanescent field corresponding to 
the nonhcrmitic part of our Morse problem. This is experimentally only a challenging 
possibility for the time being. 



References 

[1] H.C. Rosu, O. Cornejo-Perez, R. Lopez-Sandoval, "Classical harmonic oscilla- 
tor with Dirac-like parameters and possible applications", J. Phys. A 37, 11699 
(2004). 

[2] H.C. Rosu and R. Lopez-Sandoval, "Barotropic FRW cosmologies with a Dirac- 
like parameter", Mod. Phys. Lett. A 19, 1529 (2004). 

[3] F. Cooper, A. Khare A, R. Musto, A.Wipf, "Supersymmetry and the Dirac equa- 
tion", Ann. Phys. 187, 1 (1988). See also, C.V. Sukumar, "Susy and the Dirac 
equation for a central Coulomb field", J. Phys. A 18, L697 (1985); R.J. Hughes, 
V.A. Kostelecky, M.M. Nieto, "Susy quantum mechanics in a first-order Dirac 
equation", Phys. Rev. D 34, 1100 (1986); Y. Nogami and F.M. Toyama, "Super- 
symmetry aspects of the Dirac equation in one dimension with a Lorentz scalar 
potential", Phys. Rev. A 47, 1708 (1993); M. Bellini, R.R. Deza and R. Mon- 
temayor, "Mapeo en una mecanica cuantica SUSY para la ecuacion de Dirac en 
D=l+1", Rev. Mex. Fis. 42, 209 (1996). 

[4] E. Witten, "Dynamical breaking of supersymmetry", Nucl. Phys. B 188, 513 
(1981). 

[5] G. Leveque et al, "Atomic diffraction from nanostructured optical potentials", 
Phys. Rev. A 65, 053615 (2002). 
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Figure 1: The real part of the bosonic mode W2,nonherm(y) f° r x £ [0, 3] and the 
nonhermiticity parameter in the range K £ [0,2]. The rest of the parameters 
have the following values: A = 1, B = 2, a = 0.5, K' — 2, and a\ t 2 = 1 m all 
the figures of this paper. 



Figure 2: The imaginary part of the bosonic mode W2,nonherm(y) for the same 
range of x and K. 
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Figure 3: The real part of the fermionic mode ^fi, n onherm{y) f° r x £ [0, 3] and 
the nonhcrmiticity parameter in the range K G [0, 2]. 



Figure 4: The imaginary part of the fermionic mode wi. non /j erm (y) for x G [0, 3] 
and KG [0,2]. 



